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Abstract 

We propose a microscopic theory of interaction of long wave molecular 
phonons with electrons in fullerides in the presence of disorder. Phonon re- 
laxation rate and frequency renormalization are discussed. Finite electronic 
bandwidth reduces phonon relaxation rate at q = 0. Electron-phonon cou- 
pling constants with molecular modes in fullerides are estimated. The results 
are in good agreement with photoemission experiments. 
PACS numbers 79.60.Bm, 63.20.-e, 71.38.+i 

I. INTRODUCTION. 

Superconducting fullerides are a new type of materials, where the electronic bandwidth is of 
the same order as the frequencies of intramolecular modes. . Nonadiabaticity of electrons 
measured by the ratio of characteristic phonon frequency uj to the Fermi energy Ep is not 
small. The phonon frequencies are high u < 0.2 eV and bare Fermi energy is low Ep < 0.2 
eV i. 

In the past years several different calculations of the electron-phonon coupling constants 
have been reported for fullerides 0-0. Some of them yield the strongest coupling with the 
high frequency Hg modes with a moderate electron-phonon coupling, A < 0.5. on the other 
hand picket et al predicted the strongest coupling with the high frequency Ag(2) mode and 
A ~ 3. Similar conclusion has been reached in Ref. [§. The difference in calculated coupling 
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constants is quite remarkable, and may result in a qualitatively different understanding of 
the nature of superconductivity of fuUerides. Therefore, the experimental determination of 
A is required P-[T2[]. 

Recently Raman spectra for metallic fuUerides at low temperature has been reported in 
Ref. 0,0. The linewidth have been analyzed using Allen's formula for the decay rate of 
the phonon into electron-hole pair 7 averaged over all phonon momenta [ll3| : 



7 ^ (1) 

A^(0) is the density of states on the Fermi level, A is electron-phonon coupling constant, k 
is degeneracy of the phonon mode. 

It is well known that the phonon lifetime is determined by the parameter qvp/oj, vp is 
Fermi velocity, and uj is the frequency of optical phonon. It means that Allen's formula 
for the phonon linewidth does not work for optical phonons in g limit. In the clean 
single-band system it is not possible for any optical g ^ phonon to decay into electronic 
excitations because of the conservation of the momentum and the energy. 53n(g, cu) = for 
qvp <^ This result is based on the Ward identity and is independent on vertex 

corrections (Eqs. (4.20) and (4.21) of the Ref. |T^). Moreover high frequency phonon can 
not decay into electron-hole pairs if u; > w, w is electronic half bandwidth. 

It should be pointed out some special comments concerning the pentagonal pinch Ag{2) 
mode. It shows only a little broadening with doping. The authors in Refs. ||9|-[TI| conclude 
that it is the manifestation of the weak coupling with this mode. As it has been mentioned 
by Gelfand a g = mode shifts all the energy levels on all molecules in the solid 
by the same amount, and therefore leads to only diagonal elements between band states for 
the deformation potential. The q = Q Ag modes are thus uncapable of decaying into an 
electron-hole pair, no matter how strong the electron-phonon coupling is. 

Thus, because Fermi velocity vp is small and the frequency of intramolecular modes are 
high, the ratio vp/uj is small and formula for the phonon lifetime in adiabatic limit — 
does not work. The finite contribution to the phonon lifetime for c<j — appears due to 



impurity scattering and orientational disorder and the violation of the conservation of 
the momentum. 

In this paper we analyze the phonon relaxation rate and the renormalization of the 
phonon frequency for g = due to electron-phonon interaction in the presence of disorder 
and taking into account finite electronic bandwidth. We take into account the effect of 
disorder in terms of relaxation time r and adopt Fermi-liquid description EpT ^ 1. 

II. EFFECT OF DISORDER. 



We describe electron-phonon interaction in fullerides by the standard hamiltonian |T^. It 
describes the interaction of tiu electrons with Ag and Hg intramolecular modes: 

k,cr,i k,q,a,i <? 

where the first term is the kinetic energy of the electrons in threefold degenerate ti„ band, 
c\^^ is the creation operator of the electron with momentum fc, spin a and orbital index 
i{i = 1,2,3), 6j is the creation operator of the phonon with momentum q. Here we take 
into account momentum dependence of the coupling constant explicitly. For intramolecular 
modes this dependence is weak, but as we discuss later this dependence is responsible for the 
finite contribution of the electron-phonon coupling to the phonon relaxation rate at g = 0. 

Note that the fine structure of Hg phonons is usually neglected for the analysis of the 
relaxation rate with Allen's formula [1^. We also neglect strong degeneracy of Hg modes. 
This assumption is quite reasonable if electronic relaxation time is large ur >> 1 and if 
crystal field effects are strong and split of the fivefold degenerate modes is strong [1^. It is 



clear because the nodiagonal elements of the electronic Green's function appear only due to 
impurity scaterring and are small if ujt >> 1. It is important, that the interaction constant 
with single Hg submode is strongly momentum dependent. 

Phonon relaxation rate and frequency renormalization are determined by the real and 
imaginary parts of the polarization : 
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U{q = 0,uj) = i [ T{k, J + J - uj/2)gkG{k, J + uj/2)G{k, J - (3) 
J (27r) 

The equation for the vertex has the form (Fig.l): 

r(fc, J + J - uj/2) = gk + ni^/{2Tif [ \u{p - k)\^Gip, J + cu/2) (4) 



G(p,cj -cu/2)r(p,cj +cj/2,cj -oj/2)dy 

wheie u{p — k) is the potential of the single impurity, G{k,uj) = 1/{uj — ^ — Ti{uj)) is electronic 
Green function, averaged over impurity |]TB[, ^{uj) ~ ""^2}^' ^ electronic relaxation time, 



riim is concentration of impurities. 
We define the function: 

P{k, J + cu/2, tu' - Lu/2) = r(A;, + w/2, tu' - lj/2)G(/c, + ^l2)G{k, J - uj/2). (5) 

This function satisfies the equation: 

P{k, J + cj/2, J - cu/2) = ^(A;, J + cj/2)G(A;, cu' - ujl2){g^ + (6) 
ni^/(27r)=^ / - A;) |^P(p, cj' + cu/2, cu' - uj/2)(fp) 



Main contribution to the integrals appears from the momenta near the Fermi surface k ^ kp 
and we can expand and \u{k — p\'^ in spherical harmonics on the Fermi surface [jl9 



9k = Y.9L<pL{k) (7) 
Kp-A:)P=E<^^WrL,L'0L'(p)* (8) 

L,L' 

For the sake of simplicity we suppose that ^' = 5^ l'^l- The equations for the relaxation 
times have the form 1/r = 27rA^(0)?2jmro, I/t/, = 2TiN{Q)nimX' where A^(0) is the density 
of state on the Fermi level. Note that gL=o ^ S'l^^o for niodes. On the other hand for 
fivefold degenerate Hg modes we expect strong k dependence of the coupling constant. 
We define the set of functions IS.i{J ,uj): 

J2gLMk)M^' = nim/{2nf J \u{k ~ p')\^P{p' ,uj' + u/2,u' - uj/2)d^p, (9) 



and derive the equation for Kl{uj\uj): 

Y,9l(1)l{1)^l{uj\uo) = niml{2'KfY.^l + Km{uj\uj) I \u{k-p)\^gM(l)M{p) (10) 

L M 

G{p, J + cj/2)G'(p, J - u/2)d^p. 
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Integrating out the angles in Eq.(lO) and taking into account Eq.(8) we obtain: 

Kl{'^\uj)= i/tl — , , Icu'l < (11) 

UJ + 1/tI 



where l/r£ = 1/r - 1/tl 

Note that for L = Ao(ci;', cj) = i/ruj. The largest term in the expansion of the coupling 
constant go does not contribute to the g = phonon relaxation rate. Substituting Eq.(5) to 
Eq.(3) and taking into account Eqs.(9),(6) and (11) we obtain: 

n(0, uj) = -2i ^ , • (12) 

L^O ^ + V 

Here we take into account that / dcu' {Ti(iu + uj' /2) — Ti(iu — uj' /2)) = 0. 
As a result we obtain the formula for the phonon relaxation rate '^{uj): 

7(^) = -5>n(0, uj)=2}_^ (13) 

It follows from the Eq.(13) phonon relaxation rate at g — * is determined by the parameter < 
gl > — < Qk where < .. > is average over Fermi surface. This formula is strongly different 



from Allen's formula |13|. (i) Phonon relaxation rate is proportional to the averaged over 
Fermi surface /c-dependent component of the electron-phonon coupling constant. Phonon 
relaxation rate due to electron-phonon coupling is equal to zero if coupling constant is 
independent of the electronic momentum k. {ii) Phonon relaxation rate is proportional to 
the impurity scattering relaxation rate of electrons at low temperatures 1/r*. Therefore, 
momentum dependence of the electron-phonon interaction is responsible for the finite Raman 
linewidth. 
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It should be pointed out, that similar formula for the relaxation rate of the optical 



phonons in metals was derived from kinetic equation in Ref. and Green's function 
technique Note that formula (13) is different from that derived in Ref. [^]. New term 
proportional to A(uj',iu) appears in the equation for n(co') due to correct average of the 
vertex over impurities. Neglecting this term one can derive the same formula for relaxation 
rate as Eq.(18) of Ref. 



Extensive numerical calculations of the phonon lifetime, using spherically symmetrical 
coupling have been performed in Ref. [^. It has been shown that diagonal component 
of the polarization is site dependent in disordered phase. This fact is in agreement with 
formula (13). Because Hg modes are not spherically symmetrical the interaction with the 
five split submodes will have large L ^ harmonics on the Fermi surface even in the case 
of spherically symmetrical bare interaction. 

III. BANDWIDTH EFFECT. 

In superconducting fullerides there are a number of molecular modes with the frequencies of 
the order of bare bandwidths. These are pentagonal pinch mode Ag{2) uj ~ 1500cm~^ and 
four Hg modes with uo ~ 1200 — 1600cm~^. Because of conservation of energy these modes 
cannot decay into electron-hole pair in the clean system. Note that in the limit of w <^ u; 
phonon relaxation rate is equal to in the lowest order in coupling constant. 

We use Eqs. (3) and (4) for the polarization and lorenzian form of the density of states 
to take into account the finite bandwidth: 

where w is effective half bandwidth, u is orbital degeneracy. For the ti„ band z/ = 3. Using 
Eq.(14) we can derive the equation for electronic self-energy averaged over impurities in 



ladder approximation |23|: 



where x = uToriim/w'^ = 1/2tw is dimensionless concentration of impurities. 

Integrating out the angle in Eq.(lO) and taking into account Eq.(8) we obtain the formula 
for Al{uj' , uj): 

where xl = uV Lriim/w'^ = 1/2tlw. We have used here integral equation for the electronic 



self-energy in ladder approximation [^,^. Note, that Eq.(16) is equivalent to the Eq.(ll) 
if = — 2|^' Equation for the L component of the polarization has the form: 

-2tglu f ^{y-u;/2)-Eiy + u;/2) 

"^^""^ - J ^^.;+^(S(y-^/2)-S(y + ../2) ^^^^ 

Taking into account that EpT ^ wr ^ 1 we obtain: 

= XW^ ^—TT^ (18) 

^ ^ UJ + tWUj/\uj\ ^ ' 

Substituting Eq.(18) into Eq.(17) and integrating out y we derive the formulae for imaginary 
and real parts of the polarization: 

5^n(-) = E ^/if^ W 1^ (1 + ^-If] - -'^'-^ + 1) (19) 

7(07) = -53n(^) = E ^^W^T%T7V^^^^ + ^^/^)') + ^ ^'"^^^ (^/«^)/«^) (20) 

Eq.(20) reduces to Eq.(13) in the large bandwidth limit iojw <^ 1. In the opposite limit 
uj/w ^ 1 the relaxation rate is strongly reduced: 

.M^-anM^E^^^t^ (21) 



IV. CONCLUSION. 



In conclusion we analyze the experimental data on the Raman scattering in fullreides PJTO[ 
using the correct formula for g = phonon relaxation rate. Unfortunately, direct estimate of 
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the coupling constant is practically impossible. It requires exact form of angular dependence 
of electron-plionon coupling constant on the Fermi surface and electronic relaxation rate 1/r. 
However, if we assume that g1 ~< g"^ >~ Aco'/A^(0) and use the value for Hg{l) mode from 



the photoemission experiments 1T1||T^ , we can calculate coupling constants for another 7 Hg 
modes using the formula: 

7^hJ ~ A,/A, (22) 

If we suppose that Ai/A^(0) ~ 0.02eV |ini[l2| for Hg{l) mode, we obtain the coupling 
constants Aj/A^(0) for other 7 Hg modes (Table 1). Note that Eq.(22) is valid only for Hg 
modes, because angular dependence of the coupling constant on the Fermi surface for Ag 
modes is strongly different from that for Hg modes and we do not expect the cancellation 
of angular factor in Eq.(13). From the Table 1 we can conclude: 

• Using Eq.(22) and the experimental Raman linewidths we obtain coupling constants 
for Hg modes. They are in good agreement with photoemission data. Note that Allen's 
formula underestimates the coupling constants by the order of magnitude for the most 
of the Hg modes. 

• The difference in coupling constants for Hg{2) and Hg{3) modes is probably connected 
with the fact that in the analysis of photoemission spectra of Cqq the interaction with 
Ag{l) mode has been neglected |[TT|JT^ . 

• The difference in estimated constants for Hg{7,8) modes is due to frequency depen- 
dence of electronic relaxation time r. Because the interaction with low frequency 
modes is quite strong we expect strong frequency dependence of r and Eq.(22) is not 
valid. 

• Due to high symmetry of ^4^(1, 2) modes angular dependence of the coupling constants 
is weak and Eq.(22) does not work. 
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It should be pointed out that frequency renormahzation of these modes is not due to 
effects considered in the paper. Indeed, the downshift of Ag{2) mode is about 6 cm~^ per 
elementary charge on the Ceo- If "we suppose that this downshift is due to interaction of 
phonons with band electrons one should expect the maximum of downshift near the half- 
filled band (x = 3) and the absence of the downshift for x = 6. In an isolated molecule there 
is also a frequency renormahzation when molecule becomes charged. Theoretical estimates 
of the frequency shift due to the charging of Ceo molecule are in a reasonable agreement 



with experiments ||24 



We have estimated coupling constants of the conducting electrons with the molecular 
phonons in superconducting fuUerides from Raman experiments. The results are in good 
agreement with that obtained from photoemission measurements. Note that these constants 
with proper account of polaron effect lead to correct values of Tc, isotope effect and pressure 
dependence of Tc{P) [0. 

We highly appreciate enlightening discussions with A.S. Alexandrov, N.M. Plakida, J. 
Annett, E.G. Maksimov, O.V. Dolgov and D. Mihailovic. One of us (V.V.K.) thanks Russian 
Foundation for Basic Research (Grant 97-02-16705) and Slovenian Ministry of Research and 
Technology for financial support and Dragan Mihailovic for hospitality. 
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TABLES 

TABLE L Coupling constants obtained from Raman measurements using Allen's formula (AF), 
Eq.(13) and from photoemission experiments (PES). 





UJ 

(cm^^) 


7 m 

(cm^-*^) 


A/iV(0)(eV) 
AF@ 


A/iV(0)(eV) 
Eq.(13) 


A/iV(0)(eV) 
PES 


A/iV(0)(eV) 

PES m 


H{1) 


270 


20 


0.048 


0.020 


0.019 


0.020 


H{2) 


432 


21 


0.020 


0.021 


0.040 


0.038 


H{3) 


709 


8 


0.002 


0.008 


0.013 


0.019 


H{A) 


773 


10 


0.003 


0.010 


0.018 


0.018 


H{5) 


1100 


11 


0.001 


0.011 


0.012 


0.009 


H{Q) 


1248 


10 


0.001 


0.010 


0.005 


0.001 


H{7) 


1425 


46 


0.004 


0.046 


0.017 


0.000 


H{8) 


1572 


42 


0.003 


0.042 


0.023 


0.000 
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FIGURES 

FIG. 1. Equation for the vertex function r{k,uj' + uj/2,u;' — uj/2) 
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FIGURES 




FIG. 1. V.F. Aksenov, V.V. Kabanov 
Electron-Phonon Interaction and Raman Finewidth in Superconducting Fullerides. 
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